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Abstract: Relativistic axions are good candidates for the dark radiation for which there 

^ \ are mounting observational hints. The primordial decays of heavy fields produce axions 

^^ \ which are ultra-energetic compared to thermalised matter and inelastic axion-matter scat- 

00 \ tering can occur with EcoM ^ T^-, thus accessing many interesting processes which are 

otherwise kinematically forbidden in standard cosmology. Axion-photon scattering into 
^^ ■ quarks and leptons during BBN affects the light element abundances, and bounds on over- 

fT^ ■ production of ^He constrain a combination of the axion decay constant and the reheating 

temperature. For supersymmetric models, axion scattering into visible sector superpart- 

ners can give direct non-thermal production of dark matter at Ty ^ Tjreezeout- Most 
^ I axions — or any other dark radiation candidate from modulus decay — still linger today 

» I as a Cosmic Axion Background with Eaxion ~ O(IOO) eV, and a fiux of ~ 10^ cm~^s^^. 
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1. Introduction 

Conventional cosmology begins with the Hot Big Bang. During this period, the energy 
density of the universe lies entirely in thermalised Standard Model degrees of freedom, and 
the only particles with energies E ^ T are non-relativistic. 

However there are mounting hints that this standard picture of the early universe is 
not fully correct. In particular, there may exist dark radiation, i.e. non-Standard Model 
particles that are relativistic at the times of Big Bang Nucleosynthesis (BBN) and at the 
formation of the Cosmic Microwave Background (CMB). Including the HST measurement 
of the Hubble constant [1], the latest observational values from WMAP, ACT, SPT and 
Planck are N^jf = 3.84 ±0.40 (WMAP9, [2]), 3.71 ±0.35 (SPT, [3]), 3.50 ±0.42 (ACT, [4]), 
and 3.62 ±0.25 (Planck, [5]).^ While by no means conclusive, these hint at a non-Standard 
Model value of -/Vg// > 3.046, and the consequent existence of dark radiation. 



^Without including direct measurements of Ho, the determinations using only CMB and BAO data are 
N^ff = 3.55 ± 0.60 (WMAP9), 2.87 ± 0.60 (ACT), 3.50 ± 0.47 (SPT) and 3.30 ± 0.27 (Planck). Note that 
there is some overlap in the data sets in all cases and so the values should not be taken as independent. 
Also note that the straight ACDM fit to Planck gives an inferred value for Ho of (67.3 ± 1.2) km s~^ Mpc~^ 
which is low by 2 -h 3cr compared to recent direct measurements [1,6,7]. 



Dark radiation must involve a particle that is still relativistic at CMB formation at 
t ~ 3 X 10^ years, which suggests that it is extremely light. One outstanding example of 
naturally light particles are axions — either the QCD axion or more generally an axion- 
like particle. Relativistic axions can also easily be produced non-thermally through the 
primordial decays of heavy moduli or saxions. 

There is a large literature on the cosmology of axions and a recent general review 
is [8]. There is also a growing literature on dark radiation, for example [9-28]. Within 
these, [22] gives a general review of dark radiation from particle decay, and string models 
of axionic dark radiation directly relevant to this work are [20,21] (also see [17]). The role 
of relativistic thermal axions in dark radiation has been considered in [10,19,25]. 

However what this literature all tends to assume is that, once established as a dark 
radiation candidate, axions become non-dynamical. The justification for this is that axion 
scattering is suppressed by fa^i giving T/H <^ 1 and an axion population which is far from 
thermal equilibrium. In particular, the effects of their scattering off the ambient thermal 
plasma are neglected. 

In this paper we show that there is much interesting physics from the interactions of 
dark radiation axions with the thermal plasma. In particular, it allows for a violation of 
some familiar properties of early universe cosmology. For example, in Standard Cosmology, 
there are no relativistic particles with energies E ^ T since the energy spectrum of particles 
is set by the Boltzmann distribution. The ultimate abundance of a species, Y = n/ s, is 
determined either via the freeze out of stable states (as in WIMP dark matter), or by a 
conserved quantum number (as in the case of baryon number, where Yg is far in excess 
of the value given by thermal decoupling). Furthermore, any processes with a center of 
mass energy E ^ T are inaccessible: the thermal bath is simply incapable of generating 
collisions with such a high centre of mass energy and they are suppressed by e~ ' <^ 1. 

With axionic dark radiation, these properties no longer need hold. Highly energetic 
axions may arise from very simple scenarios of the early universe in which long-lived heavy 
particles of mass m$ reheat the Standard Model degrees of freedom to T^eheat ^ "t-* ) while 
also decaying to dark radiation axions with Eaxion = rn^/2. The subsequent interaction 
of these axions with the thermalised plasma allows processes at EcoM ~ \/Treheatf^^ ^ 

-'- reheat- 

This scenario is very well motivated: the existence of at least one axion (the QCD 
axion) has strong support from the non-observation of a neutron electric dipole moment, 
and string theory models in general can have many 'axion-like particles' (hence, axions) 
which may not couple to QCD but still enjoy an approximate Peccei-Quinn shift symmetry. 
This symmetry forbids perturbative mass terms for the axion, making it naturally light. 
Furthermore axions are very weakly coupled to the visible sector — for the QCD axion, 
experimental bounds require 10^ < /pg/GeV — and require large reheating temperatures 
to thermalise. For example, with fpq = lO^^GeV, the QCD axion thermalises at T > 
2 X 10^ GeV [19]. 



The existence of a long-lived massive particle species <^ which reheats the universe 
through its decay is similarly well- motivated. The relative scaling of matter [p ~ a~^) 
and radiation (p ^^ a~^) implies that, after inflation, the energy density of the universe 
tends to become dominated by non-relativistic matter particles with the longest lifetimes. 
As lifetime is inversely related to interaction strength, these are naturally particles with 
interactions suppressed by Mp^ — in the context of string compactifications these are the 
ubiquitous moduli. Particles with Planck-suppressed couplings have a characteristic decay 
rate and reheat temperature 



SttM' 



-^. TsM,reheat ~ "^ = 0-6 GcV (^^60^ j 



«' M^ 



and will come to dominate the energy density of the universe and drive reheating. The 
existence of such particles is generic in string theory or any extra dimensional theory, where 
they arise from the higher-dimensional modes of the graviton. 

However, gravity induces very general couplings and there is no a priori reason for <I> 
to couple only to the Standard Model. On the contrary, in general the field $ may also 
decay into any light axions that exist, and such axions contribute to the 'dark radiation' 
of the universe. From the decay <1> — )■ aa, the initial energy of such axions is 
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Explicit models for such decays are discussed in [9,15,20,21]. As axions are very weakly 
coupled, they do not thermalise after their original production, and their original energy 
is only diluted through the redshift due to the expansion of the universe. In particular, 
the axions remain highly energetic compared to the ambient photon plasma, by a factor 

Although we have considered the particular case of Planck-coupled moduli above, the 
same picture will occur whenever reheating is driven by the decays of long-lived massive 
states (provided they have a decay mode to axions). Other examples could be topological 
field configurations, string winding states, Q-balls etc. The key feature is simply that 
m$ > Treheat ~ [H'^MpiY/* ~ (P^Mp^)^/"^. In this case the energy of hidden axions 
produced via direct 2-body decays of <!> is again very much greater than the particle energies 
in the Standard Model thermal plasma. 

We also note here, that while we focus on axions in this paper, many of the considera- 
tions also apply to other dark radiation candidates such as hidden photons. The important 
point is only that the particles are produced from modulus decay, giving them a highly 
energetic spectrum compared to that of the thermalised plasma. 

The fact that the great majority of axions never interact with the ambient photon 
plasma does not mean no interactions occur. Although such interactions may be rare, 
they are important because the centre of mass energy of the interaction is far higher than 
can be achieved by purely thermal processes. Axion-photon scattering is then the only 



way for some interactions to occur, as they are kinematically inaccessible within standard 
cosmology. 

In this paper we study the effects of the scattering of highly energetic dark radiation 
axions off the ambient Standard Model plasma. We focus our discussion on two processes: 

1. Inelastic scattering of axions off thermal photons into Standard Model degrees of 
freedom: o + 7 — )• qiQi, where qi can be either a quark or a lepton. These processes 
can be important during and after Big Bang nucleosynthesis, and may be constrained 
by their effects on the light element abundances. 

2. The non-thermal production of dark matter, after the conventional freeze-out tem- 
perature, by the scattering of axions off the thermal plasma and into sparticles, 

a + i/g^qq*- 

The paper is structured as follows. Section ^ reviews the general formalism of moduli 
decay to give both reheating and dark radiation. Section ^ discusses axion scattering and 
how it can give rise to energy injections during BBN. Section describes the generation 
of dark matter from scattering to supersymmetric particles, and in section ||, we estimate 
the present day relic flux of relativistic axions. The conclusions summarise and describe 
possible extensions of this work. 

2. Generalities of Moduli Decay 

Throughout this paper, we assume that post-inflationary reheating is driven by the decay 
of a heavy modulus of mass m^ and decay rate T. For a Planck-coupled modulus, T and 
771$ are related by 

where k is an 0(1) constant and Mpi denotes the reduced Planck mass. The modulus 

lifetime is then 

_ 47r Mil 



K^ rwi 



(2.2) 
''(J) 

Such a modulus is long-lived and we assume that prior to its decay it dominates the energy 
density of the universe. In the bulk of this paper we will use the simultaneous decay 
approximation in which all moduli are assumed to decay precisely at time r, although for 
determining the present-day axion energy spectrum in section |5| we will be more precise. 

Prior to decay, the universe is matter dominated and a ~ t^'^ . As H decay = §1 during 
matter domination, we have 

2 3 

We assume the modulus decays with a fraction Ba to an axion (constituting dark radiation) 
and a fraction 1 — Ba to the visible Standard Model sector. We will assume that the dark 



radiation component is single species and thus only consists of this axion. The initial axion 
energy density is then 

PO,axion = BaX SH^^^^yMpi , (2.4) 

and the initial Standard Model energy density is given by, 

Po,SM = (1 - Ba) X 3Hl^yMl,i . (2.5) 

As psM = ^g*{T)T^, with g*{T) the effective number of species, we have 

TreHeat = ( ^^ .(1 " i?a)i^La,M|,,) ^^^ = K f .^"/"^ , ) ^^^ ^ • (2-6) 

The initial axion number density at time of reheating is 

'^^''''"'^'' - ^a X — - -^^^ ■ (2.7) 

In the non- interacting limit, the energy density of the relativistic axions evolves with the 
size of the universe as 

Paxiony-ti) — PO,axion I ~75" / ' \^'^) 

whereas the energy density of the Standard Model plasma evolves as 

PSM(R) = P«.SM{S^f{^)\ (2,9) 

such that the comoving entropy S ~ g^fT^B? is conserved during the expansion of the 
universe. As a result the fraction of energy density in Standard Model degrees of freedom 
gradually increases as the universe expands and (?* decreases. The hidden sector branching 
ratio Ba relates to the experimental observable AN^jf = ^eff — 3.046 as (e.g. see [20]), 

A AT "^ <* I g*{-'^i/ decoupling) \ /„ i r\\ 

ANeff = —- — [ , ^ . (2.10) 

t i^ — Oa \ g*[-Lreheat) J 

After modulus decay, the Standard Model sector rapidly attains thermal equilibrium 
and the number densities and energy distribution of particles in it are set by the Boltzmann 
distribution. The axion decay mode ^ ^ aa however produces monoenergetic particles 
with an energy Ea = m<j>/2. This energy redshifts as R^^ but does not otherwise change.^ 

As a numerical example we will use throughout this paper, if m^ = 5 x 10^ GeV and n = 
1, then a hidden sector branching ratio of Ba = 0.145 corresponds to a reheat temperature 
of Treheat = 1-00 GeV. This branching ratio is chosen to correspond to the Planck + HST 
central value of Neff = 3.62, using g^{Treheat) = 61.75 and g^:{Ti, decoupling) = 10.75. The 
choice of a decaying modulus mass in the region m$ = 10^ -^ 10 GeV is motivated by the 
appearance of this scale in various string scenarios [29-31]. 
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Figure 1: The three distinct energy scales and their evolution with time: the energy of the 
relativistic axions, the energy of the thermal Standard Model plasma, and the centre of mass energy 
for scattering between the axions and the thermal plasma. These can be separated by several orders 
of magnitude. 



In this paper, our main interest will lie in the scattering of the relativistic axions off 
the thermal Standard Model plasma. In figure |l] we illustrate the approximate distribution 
of energies as the universe evolves, and how they differ from standard cosmology. 



3. Axion-Photon Scattering Constraints from BBN 



In this section we compute the inelastic scattering rate for an axion off thermal photons and 
determine the resulting constraints due to energy injections during BBN. For concreteness 
we consider the process, 



a + 7 — )• e + e"^ , 



(3.1) 



as an example of the general class of processes including axion scattering into quarks, 
a + 7 — )■ gj + ^i, or any lepton a + 7 — )■ /~ + Z^. The generalisation to different final states 
is trivial. 



^The monoenergetic spectrum is a feature of the instantaneous decay approximation. In practice the 
decays will occur at different times and the spectrum will be smeared by redshifting. 





Figure 2: Feynman diagrams contributing to a + 7 — > e e+. The third diagram involves an 
additional factor of aEM and so we neglect it. 

3.1 Calculation 

The potential Lagrangian interaction terms are 

7H — 

Axion-electron-electron: Ce—T^aip^xp , 

J a 

Photon-electron-electron: eA^ipj^^ , 

Axion-photon-photon: c^-^—e^i,xpF^^F^ = ga^^aEi-'R, (3-2) 

where Ce and c^ are model-dependent constants. We shall henceforth take Ce = 1. The 
process (^]^) proceeds via the Feynman diagrams shown in figure Q. As the third diagram 
involves the 077 vertex which is suppressed by an additional factor of ^^, we shall not 
consider it further here.^ 

We aim to compute the axion scattering rate, 

r„^^e+e- = r = {nav) , (3.3) 

by first computing av, averaged over final state momenta and initial directions, for an 
axion of energy Ea scattering on a photon of energy E^, and then subsequently performing 
the thermal averaging over E^. 

The axion and photon 4-momenta are 

kphoton = iE^,E^cos0,E^sm9,O). 
where we have assumed rua ^ Ea- The centre of mass energy is then given by E^^j^ = 

'^kaxion ■ kphoton = 2-Ea£'^(l - x) with X = COS 9. 

Denoting the amplitudes for the first two diagrams in figure |2| by Ma and Mb respec- 
tively, we have 

I . . . . i9 2e^m^ fu t \ ,„ ,, 

|-Ma + -MsP = ^^(- + --1J, (3.4) 



^The presence of oem can be understood by starting with the 'natural' form of the gauge kinetic action, 
-r4 — FfijjF'^'^ + -j^Ffj^vF'^" , and then canonically normalising the kinetic terms. 



where the first term arises from |A^/ip, the second from |A^Bp, and the third from the 
cross term. Here t = ki ■ k^ and u = ki ■ k^. Integrating over the e^e~ phase space, we get 

1 e'ml (^ ( e'+p \ ^p\ 

^-4Ki?,K-^,|2vr/2 y''\E'-p') 2E')- ^'-'^ 

Here E = EcoM /'^ is the energy of the electron in the centre of mass frame and p = 
y/{E')'^ — TTig is the momentum of the electron in the centre of mass frame. 

We can average av over the initial photon direction by performing the x-integral, 
^ J_i dx, where A = ^"^ accounts for the threshold energy for e~^e~ pair production. 
The relevant integral is 

(av) = j07fl\^dx In (^T^x + VI - x - X) - In f^T^x - Vl-x-X) 

-2 1-x ' ('•') 

which can be performed to obtain 

M = T^ -4^1nVA-^^/2^A-(2-A)ln(^/2-^/2^A 
167r/2 \^ 2 ^ ^ ^ V 

+ f2 + ^Vnf^/2 + ^/2^ 

Finally, we need to integrate {av) over the photon spectrum using the thermal number 
density, 



dn 



g E^dE^ 
2^~^ 7' 

g T — 1 



with 5 = 2. Using E.^ = -^^, we can rewrite ^^^y ^ = -^^ — W^ ' ^^^h A running from 

eBAT-l 

to 2. We then end up with, 

T. , ^ 9 "^e e^ r^ 1 1 / 5A, , 5 , ^ 

r = inav) = — TT-^^ TTT / "A 5 —7 — — m A ^v2 — A 

^ ^ 2^2 El 27: fl Jo ^ J5I _ 1 aH 4 V2 

-(2-A)ln(^/2-^/2^A)+(2 + y)ln(^/2 + V2^A)') , (3.7) 
which can be evaluated numerically as a function of the relevant parameters. 
3.2 Bounds on highly relativistic axions from Big Bang nucleosynthesis 



Equation (3/7) gives the scattering rate of energetic axions off thermal photons, and in this 
section we will show how this scattering rate — despite its relative smallness in relation to 
the Hubble parameter in the early universe — can be large enough to affect the primordial 
^He abundance. 



After the decay of the modulus ^ into Standard Model degrees of freedom and axions, 
the number density of axions is diluted with the expansion of space-time. Occasional scat- 
tering of axions off the thermal plasma into 'primary' scattering products such as e"'"e~ or 
quark/anti-quark pairs, quickly thermalize by producing a large number of electromagnetic 
(i.e. 7, e^) or hadronic (e.g. p, n) 'secondary' particles. The secondaries may in turn affect 
the light element abundances in the early universe through a number of different processes, 
as described in detail in e.g. [32-35]. This way, relativistic axions may in principle be 
constrained by Big Bang Nucleosynthesis bounds on nuclear abundances. 

3.2.1 Inelastic axion scattering and decaying neutral particles 

Unfortunately, to our knowledge there exists no dedicated study of what constraints may 
be inferred from light element yields during BBN on a highly relativistic species which 
scatters inelastically off the thermal plasma. However, there are well-known constraints 
on massive particle species, here denoted X, decaying with life-time tx during and after 
BBN. These bounds are often quoted in terms of 

ex = Mxix , (3.8) 

where ^x = nx/n^ is the relative abundance of X-particles to photons prior to the decay, 
and Mx is the mass of X. The quantity ex measures the total energy deposited per photon 
once the entire species X has decayed. Axion scattering off the thermal plasma differs from 
a decaying particle in four ways: 

1. The scattering rate, and thereby the effective 'life-time' (r = T""*^), depends on the 
temperature of the thermal photons as well as the energy of the axion, which both 
decrease with time. This gives rise to a particular time-dependence for the scattering 
rate into each primary final state. 

2. While the decay products of a decaying massive particle are traditionally assumed 
to have negligible total momentum with respect to the rest frame of the plasma, 
the primaries of axion scattering will necessarily appear boosted with respect to this 
frame. Yet, the average total momentum over many inelastic axion scatterings will 
vanish, and we expect this point to be of less importance to the applicability of the 
bounds. 

3. The role of the factor of Mx in ex which determines the total energy of the decay 
products is in the case of axion scattering played by the time-dependent axion energy 

4. Finally, while the number of X-particles is greatly reduced at times t/rx ^ 1, the 
number density of axions per comoving volume does not change by much for the 
scattering rates that we will discuss. In other words, axion scatterings are relatively 
rare events, and most relativistic axions pass through BBN unperturbed. 



Despite these obvious differences, we find that at the level of the energy deposition 
history — which is in the end what is constrained by the BBN analysis — each axion 
scattering channel can be well approximated by a species of decaying particle. Thus, the 
effects of the axion during BBN can be modelled by a collection of decaying particle species, 
each decaying only to a particular final state qiQi with life-time Tj and deposited energy (per 
photon) ej. The corresponding hadronic branching ratio for each of these scattering modes 
will be either 1 (for quark and heavy lepton primaries) or ~ 10~^ (for electron primaries). 

In this section, we impose known bounds on the energy deposition of each individual 
such species during BBN, and interpret the bound as a constraint on the axion decay 
constant and modulus mass m$. These bounds are necessarily conservative as the effect of 
each possible primary is considered separately, yet we will see how they can give important 
constraints on the axion decay constant. 

Specifically, constraints on a neutral decaying particle during (and after) BBN are 
conveniently phrased in terms of bounds on 

ex = - dt'Mxix{t')= dt'Mx\Tx\U{t'). (3.9) 

JO Jo 

given the life-time tx- For axion scattering into a particular final state, here denoted qq, 
the corresponding bounds can be phrased in terms of 

oo 

dt'Ea{t')\Ta^,^qq{t')\Ut'). (3-10) 

to 

where ^a = na/n^ denotes the relative number density of axions to photons. Here io 
denotes an initial time before which no bounds can be placed on the energy deposition, 
and during radiation domination, Ea{t') = Eq [-^) ■ Thus, the deposited energy from 
axion scattering into qq-paiis from the time fg up to t is given by. 



6,(t) = Eo j\t' (^) \Ta,^,q{t'Ma(t') . (3.11) 



Over a period in which the effective number of species remains the same, the temperature 
of the photon gas depends on time as T = Trhy/trh/t, and the re- heating temperature is 
determined from the mass of the decaying modulus (and thereby Eq) as in equation (p.6|). 

Thus, the energy deposition for axion scattering exclusively into qq can be found 
explicitly, and an example thereof is shown in figure H for scattering into bottom quarks. 
From this figure, we note that the energy deposition history from axion scattering is very 
similar to that of a decaying particle species, exit) = Mx{£,x{to) — Cx{t)). This allows 
us to map each axion scattering mode into a species of neutral decaying particles with 
life-time, mass and initial number density determined by optimising the fit of the energy 
deposition profile to that of the corresponding axion scattering. 
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Energy Deposition Histories 




Figure 3: Energy deposition histories ea{t) for axion inelastic scattering to bottom quarks (solid 
lines) as compared to the best-fit profile of ex (t) for a decaying particle species (dashed lines) . For 
both solid curves, fa — 10^ GeV and ^N^jf = 0.57. In blue, m$ = 5 • 10^ GeV, and in green 
TO$ = 5 • 10^ GeV. 



3.2.2 Observational constraints from BBN 

During BBN (10^^ ^ ^ ^ 10^ seconds), the strongest constrains apply to processes which 
increase the relative number of neutrons to protons as these neutrons will in turn increase 
the Helium abundance, Yp. At the very beginning of this period, i ^ 3 • 10~^ seconds, 
even these bounds are ineffective as thermal weak interactions quickly restore the neutron 
to proton ratio. However, particles decaying later during BBN may be constrained from 
bounds on overproduction of ^He. 

The primordial Helium abundance can be inferred both from spectroscopy of extra 
galactic Hn regions by extrapolation to zero metallicity, as well as from fits of cosmological 
models of the temperature anisotropics of the CMB and the large-scale structure of the 
universe (for a recent review, see [36]). Recent spectroscopic studies have found Yl = 



BBN 



0.2528 ± 0.0028 [38], and Y 



BBN 



0.2565 ± O.OOlO(stat) ± 0.0050(syst) [37], Y^ 
0.2534 lb 0.0083 [39]. These values are compatible within la to those inferred from the 
recent Planck experiment, which estimated Y^^^ = 0.266 it 0.021 for the ACDM model 
extended with Y^^^ as a free parameter, and Y^^^ = 0.254l{];[]33 for ACDM extended 
with both Y^^^ and A^'g// (the latter parameter is then determined to N(,ff = 3.33lQ'g3). 
In this paper, we will use the constraints on decaying neutral particles from [34] which are 
based on the conservative bound of 1^ < 0.258. 



p 
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Figure 4: Conservative bounds on tlie axion decay constant as a function of tlic initial axion 
energy obtained by considering the effects of the decay channels to bh (solid) , cc (dashed) , and ss 
(dot-dashed) separately. The red, black and green curves correspond to ANeff = 0.1,0.5 and 1, 
respectively, and the areas below the curves are excluded from the constraints from [34] due to 
overproduction of *He during BBN. 

Given the resulting constraint on S,x for a given life-time, a (conservative) bound on 
the energy deposited by the scattering axion into each decay channel may be derived. The 
resulting constraints as obtained separately for the bb, cc and ss scattering channels are 
shown in figure ^ We note that these constraints probe an interesting region of the axion 
parameter space, and that dedicated studies of inelastic axion scattering during BBN may 
well be able to derive even stronger bounds on fa- 



4. Dark Matter Axiogenesis 

One popular model of dark matter is that it consists of the lightest supersymmetric particle 
(LSP) whose stability is protected by R-parity. Conventionally, the abundance of the 
LSP is determined via a thermal freeze out from the ambient plasma at a temperature of 
approximately Tfreeze-out - mLSpf^O. 

However, the scenario considered here offers a radically new production mechanism. 
Even when the thermal plasma has T < Tfreeze-out ^ fniSP-, energetic axionic dark 
radiation can directly produce supersymmetric particles through scattering off the ambient 
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Figure 5: Dark matter production from axion scattering. Here the axion scatters off a gluon, 
leading to pair production of heavy squarks which then cascade decay to produce the LSP. 

plasma. The produced susy particles will undergo a rapid cascade decay to the LSP, 
which will remain as a stable dark matter relic. A sample process leading to dark matter 
production is shown in figure 0. The rate of dark matter production is set by the axion 
scattering rate. This is an entirely non-thermal process, and can give significant LSP 
production at temperatures T <C Tfj-eeze-out- 

Let us first give the basic scaling estimates for the dark matter abundance produced (a 
proper calculation follows). For the estimate we assume the axion-plasma center of mass 
energy is above threshold for producing susy particles and neglect mass effects. We also 
use only the 'universal' 077 or agg vertex shown in figure ^, so there is no mass suppression 
of the axion-matter coupling. 

The density of plasma particles is n ~ T^ and the axion scattering cross-section is 
cr ~ ( ^-7- 1 . The total scattering rate of an axion with the ambient plasma is then 



a \2r3 
4^^ 1! 



T = {nav)^{ — ) -2- (4.1) 



As during radiation domination T r^ t ^/^, the scattering rate T r^ t ^'^ and it is clear 

that the total number of scattering events Mtotai = / Tdi is dominated by those occurring 

at early times. As the Hubble scale is -ff ~ 'm~^ ^^^ fraction of axions that interact in a 

Hubble time is 

r TMpi / a \ 2 

As this is dominated by early times, in the instantaneous-reheating approximation, we 
simply evaluate this at T ~ T^eheat to get an estimate of the total fraction of axions 
that scatter. Each decaying modulus of mass m$ generates either two axions of energies 
Ea = m<j)/2 or ~ m^/Treheat thermalised Standard Model particles. If the axion branching 
ratio is Ba, the axion abundance is 

■^^ '^ axion J^ a -'-reheat 

^ axion = '^ • 

S m$ 

As each axion scattering event to R-parity odd particles produces two dark matter particles, 
the resulting dark matter abundance is 

'dm, '^ J^ axion ^t '^ "■ i-Q \ ,1 ) " 

H rn<s>Jd ^47r/ 
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q + q 





a + q 



q + g 



mmm 





Axion 

Gluon lAAAAAIVW/ 

Gluino ,/VlAAAIWl, 

Squark ^^^^^^_ 

Quark 



Figure 6: Feynman diagrams that lead to LSP production from axion scattering off the ambient 
plasma. We have shown diagrams that involve strong interactions but there are also similar diagrams 
involving weak or electromagnetic scattering off photons or leptons. 



While this neglects factors of (7*, vr etc, it is clear that for reasonable parameter values this 
can be a phenomenologically interesting source of dark matter (note that Y^m ~ 1-5 x 10~^^ 
for a 300 GeV LSP). 

Overall, any source of axion-plasma interaction that produces susy particles can pro- 
duce dark matter. Such diagrams are shown in figure |^ We will return to the full set of 
diagrams of figure in future work. Here however we focus only on the universal diagrams 
involving s-channel photons or gluons. These are universal as the coupling agg is neces- 
sarily present for a QCD axion, whereas any interactions with quarks or squarks are more 
model dependent with unspecified coefficients. 

4.1 Scattering via Photons 



We consider the scattering process a + 7 — >• i 

axion as massless, and give the squarks a mass of rriq. In section 



via an s-channel photon. We treat the 
2| we will consider the 
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analogous process a + 5 — )■ qq* . The 077 interaction is 



•'-077 — <^5a77-t^ ■ -t* — Q 05077 C/^j/Ap-'^ -f^ • l^-*^J 



The resulting amplitude is, 

where /ci and /c2 are the incident 4-momenta of the axion and photon, and k^ and /c4 the 
4-moinenta of the outgoing squarks or sleptons. In the center of mass frame we can without 
loss of generality write 

ki = {Ecm,0,0,Ecm) , 
^2 = {Ecm,0,0,—Ecm) , 
h = iEcM,0,psme,pcos6) , 
/c4 = {EcM,0,-psm6,—pcos9), 
e = (0, cos (j), sin 4>, 0) , 

giving 

I A^ 1 2 = e^gl^y cos^ (j) sin^ 9. (4.4) 

Averaging over the initial photon polarisation then gives 



2 2 
\M\'' = -^p' sm' . (4.5) 



If the lab frame 4-momenta are 



ka = {E, 0, 0, E) , 

k^ = {E-y, 0, E^ sin 6, E^ cos 9) , 

then the resulting cross-section, integrated over final state phase space, is 



o2„2 



1 e 5a77 L 2m2 (EE^ ^ 



4EE^ 247r y EE^{1 - x) \ 2 ^ ' J ' ^ ^ 

where x = cos 9. We finally average over the relative angle between the photon and axion, 
{av) = ^ f_^ avdx. In terms of A = ^§-, this is 

2 2 
{av) = -^^^^(2{4-5X)V^^^ + 6XHn(^V2 + V2^x)-3XHnXj , (4.7) 

for the x-section (av) of an axion of energy E incident on a single photon of energy E^. 
We finally sum over the spectrum of thermal photons by integrating 

9^ f ^7 .J, _ 9^ /"' .> 1 8m 



dE^ = ^^ / d\ 



2t,2J ^E,/T_i 7 27r2 7o ^^ _ I E^ X'' ' 
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The scattering rate to a pair of charged scalars is finally 

^ 5 

2m/- 

eE\T — 1 



^ ' 2tt^ Jq ^^_^E^X^ 128 247r V ^ ^^ 

+6X^ In (V2 + ^/2^^) - 3A2 In a) . (4.8) 



We obtain the total scattering rate to sleptons and squarks by summing over all species. 
For simplicity, we assume a common mass scale for all squarks and sleptons. Writing 

6 = ly T§f^ ^6 have 

Ttotai = 6xr((7 = l) + 6x3xr((/ = 2/3)+6x3xr(g = l/3) . (4.9) 

V V V 

sleptons up squarks down squarks 

As each axion scattering event produces pairs of susy particles, each of which undergoes a 
rapid cascade decay down to the LSP, the overall LSP production rate is 

TlSP = 2 X r total . (4.10) 

4.2 Scattering via Gluons 

There is also a kinematically identical diagram where the s-channel photon is replaced by 
a gluon. This differs only by the colour factors and the conventional notation for axion- 
photon or axion-gluon couplings. The agg vertex comes from the interaction 



a a. 



■eo^^^sGfGZ'eC. (4.11) 



fpQ Wtt 
By comparison of eq. ( 4.11| ) with eq. ([4.2|), we see that we have to replace 



{9 a'Y'y) photon with | — I 



1 as\ 

strong 



and 

-ie(/c3 - A;4)^ with - igsi^ks - k4)H^ , 

„2 

where Qs is the strong coupling and a^ = |;^. Summing over the final squark colour states, 
in eq. ( ^.11| ) we replace 

2 „2 

e g^^^ with 



-^'"^^"-(7^£)| 



As all possible scalar final states are squarks, we have 

Ttotai = 12 X r, and again Tlsp = 2 x Ttotal ■ (4.12) 



16 



4.3 Dark Matter Abundances 

In the instantaneous reheating approximation, we will take TT/n^"'""'* { as a proxy for the 

^ K-'-reheat) 

total fraction of axions that scatter. The resulting dark matter number density is 

o /rr \ ■'- {^reheat) f a 1 o\ 

niSP = ^ X naxionyJ- reheat) X 77777; T • l^.idj 

^{.J' reheat) 

The axion number density is given by eq. (^J) which we repeat here, 

na.^on,0 = B. X ^H^Treheat)'^h (4.14) 

-C/a 

with Ba the modulus branching ratio to axions. 

As the entropy density is s = ^(7*T^, the dark matter abundance is given by 

^r '^LSP ^t) flaxiony J- reheat) '-{-'^reheat) /, ^ -\ 

•■^ 9*-^ ^reheat H{Treheat) 



'^{Treheat) IS Set by eq. ( |4.12D . H{Treheat) IS determined by taking 



where k is a parameter we scan over. 

We scan over values of m^ between 10^ GeV and 10 GeV, and values of k from 0.5 
to 5. The results are shown in figures |^. The figures use N^ff = 3.62, fpg = 10^ GeV and 
a common scalar mass of 1 TeV. We also set Og = 0.11 and use g^: = 61.75. It is clear that 
it is easy to obtain dark matter abundances compatible with observations. 

5. Present Day Cosmic Axion Background Flux 

As we have discussed in the previous sections, dark radiation in the form of relativistic 
axions rarely interacts with its surrounding and the vast majority of the particles pass 
through the cosmos unperturbed, leaving a present day population of relic axions which 
is only diluted by the expansion of space. In this section, we compute the flux of axionic 
dark radiation in the current era (hence, 'today'), and find that it may even dominate over 
the flux of solar axions which are currently searched for by e.g. the CAST experiment at 
CERN. We furthermore point out that the differential axion spectrum and the conversion 
probability differs between solar axions and axions arising as dark radiation from moduli 
decay. 

We re-emphasise here that while we focus for definiteness on axions, the same spectrum 
and flux will also hold for any other dark radiation candidate — for example hidden photons 
— that is produced by modulus decay. 

Let us first discuss the flux of primordial axions arising from the instantaneous decay 
approximation for the modulus $, thus constituting a mono-energetic population at all 
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Dark Matter Abundance from axion-gluon scattering Reheat Temperature (GeV) 




Mj./lO" GeV 



Mo/10^ GeV 



Figure 7: The dark matter abundance and reheating temperature as a function of k and m$ . We 
have assumed N^/f = 3.62, fpq — 10^ GeV and a common squark and slepton mass scale of 1 TeV. 



times. Immediately after the decay of <I>, the population of axions is given by equation 
(I2: 



4 ^5 



IT'axionfi 



BgK^ m% 
67r2 Ml^ 



,32 r, .a( ^^ \ 



3.7 • 10'^^ • Bate 



106 GeV/ 



cm 



(5.1) 



which we repeat here for clarity. Using the expression for T^h from equation (p.6|), the 
entropy at time of reheating is easily found to be, 



Sk* ^ I a f ^ \ ' 777 



27r^ 



45 



457r' 



3/4 _9/2 



M 



3/2 ' 



(5.2) 



PI 



where g^, = g*{Tj.eheat)- The axion abundance — both initially and at the present day 
is then 



-f axion — — ^a 



127r 1/4 V5(l-i?a)y Ml/' 



"' mf 



3.85 • 10" 



K Ba 



^1/4 (l_ 5^)3/4 V 106 GeV 



PI 
m<j) \ 1/2 



(5.3) 
(5.4) 



With the present day CMB photon density of n^ = 413 photons per cm ^ , the axion 
number density today is given by, 



Tia — ^ axiari'^ ~ ' •^^'^'y ^ axion — i.lZ * lU 



K Ba 



^1/4(1-5^)3/4 V106 GeV 
The corresponding flux of dark radiation axions is then simply given by. 

-16 



ms, \ 1/2 o , ^ 

^ cmT^ . (5.5) 



$„ = -^ = 8.39 • 10^ 



gi/4 (1-5^3/4 V106 GeV 



{- 



rrii^ \i/2 -, _2 
' s cm 



(5.6) 
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For the example we have considered throughout this paper with m$ = 5 • 10^ GeV, k = 1, 
Qi, = 61.75 and Ba = 0.145, the present day flux is, 

$a = 1.09 • 10^ s^^cm^^ . (5.7) 

This flux may be compared to that of axions created from Primakoff scattering in the 
sun, which — contrary to the isotropic dark radiation flux of equation ( |5.6|) — is always 
directional and suppressed by two powers of the axion decay constant [40], 

^ solar = {ga-ry X 10^°GeV) ^ 3.54 X 10^^cm~2g-i 

R/10i°GeV\^ ,1 , , 

= 1.91 X 10^ ( J cm^^s"^ . (5.8) 

The present day flux of dark radiation axions can easily dominate the solar axion flux for 
non-excluded values of the axion decay constant. 

Moreover, the spectrum of solar axions differ significantly from that of axions produced 
from moduli decays. The differential solar flux spectrum is given by 

-( ^'"^ 4.02xlO^°cm'V^keV- \,;,,^,,J ^ , (5.9) 



dE V vr ; /2 ■ ei?/(i.08keV)_i ' 

which peaks around 4keV. 

In the instantaneous decay approximation — which we will soon go beyond — the 
ratio of the axion energy to the CMB temperature is easily obtained. Initially, £"„ = 
m<j,/2 and T^ = Treheat- The axion energy redshifts with the expansion of the universe as 
-E-a ~ "^5 while the photon energies redshift as Ty ~ ^^3 . There are two boosts to the 
photon energy: from the reheat temperature to the time of neutrino decoupling, and from 
neutrino decoupling to now. These give 



J' / r, (T^ . .\ \i/3 /ii\i/3 



Ea 

which evaluates to 



*\-^ u decoupling) J \^J \Ea J reheat 



Ty niV^^2F ^ 1 \i/3 /^n _ R ^^l/4 /_,\i/2 



Ea 



f(l^)(^^) «•<--'- (S)' ^-^' 



Thus, in the instantaneous decay approximation, the energy of the dark radiation 
axions today is given by. 






E,{T,) = 1.97 -^ -^jj^ — , (5.12) 



where we have again abbreviated g^, = gi^iT^eheat)- In the example considered above with 
m$ = 5 • 10^ GeV, the present day energy of the axions is, 

^a (today) = 238 eV, (5.13) 
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which is an order of magnitude lower than the energy of solar axions. 

In the above discussion, we have used the instantaneous reheating approximation, 
where all moduli decayed at a time r. This produces a monoenergetic axion spectrum, 
with a single initial energy of Ea = m^/2 that is subsequently redshifted. In reality 
the moduli decay gradually, and the axion energies are redshifted by different amounts 
depending on whether they come from early-decaying or late-decaying moduli. 

We now want to determine the actual form of this spectrum. We note that this 
structure of the axion energy spectrum is essentially fixed by a time t ~ (a few)T. At 
this point all the moduli will have decayed and the universe will have transitioned to 
radiation domination. As the axions are to a good approximation non-interacting, the 
energy spectrum now can be obtained simply by redshifting this spectrum. 

To determine the structure of the energy spectrum we can calculate within a comoving 
volume: there is then no dilution in number density due to expansion of the universe. To 
convert to physical number densities, we simply multiply by a factor of R~^ - which clearly 
does not affect the form of the energy spectrum. 

For a modulus with lifetime r, the fractional decay rate is 

dJy = , givmg Jy = I\Qe ^. 

T 

The initial axion energy from a decay at time t^ is always Eq = -^, and at time t > t^ is 



Ei, — En 



2 ' 

ajtd) 
a{t) 



Axions generated between times t^ and td+dtd then lie in energies between Et and Et+dEt, 
where 

dEt = ^^Etdtd = EtH{td)dtd, 
a{td) 

and the number of axions generated between times td and td + dtd is 

dNa = ^Noe-'-^dtd = '^-^N.e-'-^-.^dEt . (5.14) 

At any time i, we then have the axion number spectrum 

dN 25 _M 1 .__. 



dEt T EtH{U 

where t^ and H{td) are implicitly related to Et by Et = Eq 



a(td) 
ait) 



The solution of eq. ( 5.15| ) requires the determination of a{t). This can be found 



numerically as the modulus gradually decays and the universe transitions from matter 
domination to radiation domination. The relevant equations are 

/9$ + 6Hpis> = , 

T 
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Figure 8: A typical axion flux per square centinietre per second, with -/Vg// = 3.62. The integrated 
flux is 0.96 X 10^ cm~^s~^. The precise location of the energy peak depends on the value of m$. 
The cutoff at low energies is a numerical artefact from the spectrum computation. 




3M2, 



Using the numerical solution of a(t), the resulting form of the axionic flux is shown in 
figure |8|. The spectrum is normalised to N^jf = 3.62. The exact location of the energy 
peak relative to the CMB temperature depends on m^, and a precise determination of 
this requires an analysis of Standard Model reheating beyond the instantaneous decay 
approximation. However for m^ ~ 10^ -^ lO^GeV it will always be more energetic than the 
CMB peak by an approximate factor of 10^. 

We can obtain analytic approximations using the simplifying assumption of either 
continual radiation domination or continual matter domination. These assumptions are 
accurate for either late modulus decays (radiation domination) or early modulus decays 
(matter domination). 

For radiation domination, as a ~ t^'^, we can relate t^ and Et by 



td 
t 



El 

Eq 



(5.16) 



giving a differential axion number density of 



dn{t,Et) 2BaNo -if|i 



dEt 



rEt 



2t 



El 

Eq 



(5.17) 
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Continual matter domination with a ~ t^' ^ gives 



dn(t,E,) 2B„JV„ -i(|i-)"'3t /E,^"^ 



This spectrum is characteristic of dark radiation from moduh decays, and may be 
searched for by experiments hke CAST that look for axion-like particles through conversion 
in a magnetic field. We note that since the Cosmic Axion Background flux arising from 
modulus decay is independent of fa, the axion-to-photon generation rate in a magnetic 
field scales as /q"^, as compared to the f~^ fall-off of the axion-to-photon generation rate 
for solar axions. 



6. Conclusions 

In this article, we have studied the rich and interesting cosmophenomenology of axionic 
dark radiation. We have focused on three disparate areas: additional energy injection 
at BBN from axion-matter scattering, the late-time production of dark matter through 
scattering of axions into susy particles, and the relic Cosmic Axion Background flux today, 
which can be comparable to or larger than the solar axion flux. Our main point is that 
there is a rich phenomenology associated to axionic dark radiation, which is by no means 
sterile. 

There are several clear directions for future work, one of which is the extension of 
our computation in section ^ to a larger number of processes for axion scattering into 
dark matter. In this initial study we principally considered the simplest case of universal 
s-channel gluon scattering. There are also many other diagrams, enumerated in figure 
|6|, which we did not consider. These diagrams may in fact be more important, as the 
077 or agg vertex is generally suppressed by an extra factor of -^ compared to axion- 
matter vertices. Axion-matter-matter couplings do have a suppression by the mass of the 
relevant matter particle, but this should be less relevant once the centre of mass energy is 
comparable to the particle masses. 

Furthermore, we have emphasised that bounds from overproduction of ^He during 
BBN place interesting constraints on the axion decay constant and the reheating tempera- 



ture. While our bounds in section 3.2 were derived for each decay channel separately and 
therefore constitute conservative bounds, a full, inclusive, analysis of the effects of dark 
radiation during BBN is likely to find even stronger constraints. 

In section 5 we pointed out that dark radiation from modulus decays has a characteris- 
tic isotropic and non-thermal spectrum at an energy scale around six orders of magnitudes 
higher than that of the CMB. While we have focused on axions, the same spectrum would 
apply to other dark radiation candidates such as hidden photons, which may also intercon- 
vert to photons in the presence of a magnetic field. If it exists, this spectrum represents a 
Cosmic Axion Background which gives an image of the universe at a time t ~ 10~^ s, and 
its significance is obvious. As the flux is not dissimilar to the solar axion flux in magnitude. 
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but isotropic in direction, it represents an exciting and worthwhile target for experiments 
such as CAST that rely on axion-photon conversion within a magnetic field. 
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